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ABSTRACT 

Aims. We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov- 
Gerasimenko and provide constraints on the mechanical properties of the cometary material (tensile, shear, and compressive 
strengths). 

Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morpholo¬ 
gies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), using two shape models computed from OSIRIS images by 
the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive 
strengths using different surface morphologies (overhangs, collapsed structures, boulders, cliffs, and Philae’s footprint) and mechani¬ 
cal considerations. 

Results. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain mor¬ 
phology: i) low-slope terrains (0-20°) are covered by a fine material and contain a few large (>10m) and isolated boulders, ii) 
intermediate-slope terrains (20-45°) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size 
boulders from <1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90°) are cliffs that expose a consolidated mate¬ 
rial and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3 -15 Pa (upper limit of 150 Pa), 

4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (up¬ 
per limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive 
strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the 
formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and 
tends to favor a formation of comets by the accrection of pebbles at low velocities. 

Key words. Comets: individual: 67P/Churyumov-Gerasimenko - Comets: general - Accretion, accretion disks - Methods: data 
analysis 


1. Introduction 

Rosetta has been orbiting comet 67P/Churyumov-Gerasimenko 
(67P) since August 2014. The OSIRIS cameras ( [Keller et al.| 
2007 j ) onboard this spacecraft have acquired hundreds of images 
of the surface with an unprecedented spatial resolution down to 
the decimeter scale ( jSierks et al.||2015| ). The images reveal a 
complex nucleus surface made of smooth and hummocky ter¬ 
rains that are partially or entirely covered by dust or expose a 
consolidated material, pits, cliffs, and fractures from the hundred 
meter scale to the decimeter scale ( [Thomas et al.|2015] ). The na¬ 
ture and origin of these terrains and morphological features are 
far from being understood, but remain of paramount importance 
to better constrain the formation and evolution scenario of the 
nucleus of 67P and comets in general. This paper focuses on the 
link between the nucleus gravitational slopes and surface mor¬ 
phology to provide constraints on the nature of the cometary ma¬ 


terial and its mechanical properties in particular (tensile, shear, 
and compressive strengths). 


Gravitational 
cometary nuclei 


slopes have only 
so far, 9P/Tempel 1 


been measured on 
([Thomas et 


al 


three 

2007[), 


8IPAVild 2 ( [Jordaetal.[2QT5l ), and 67P ( [Jordaet al.[20T5l i. While 
the slopes of 9P are between 0° and 35°, those of 8 IP and 67P 
cover a much wider range from 0° to >90°, slopes exceeding 
90° indicate overhangs. Beyond the different spatial resolution 
of the shape models used to compute the gravitational slopes for 
these three bodies, the differences between 9P on one side and 
8 IP and 67P on the other side are interpreted as an aging effect 
by [Jorda et al.[p015] ). Polio wing the scenario described by these 
authors, comets that have spent more time in the inner solar sys¬ 
tem like 9P have been smoothed and have a narrower range of 
gravitational slopes than comets that have spent less time in the 
inner solar system like 8 IP and 67P. This planation process of 


Article number, page 1 of 18 























ASlA proofs: manuscript no. 26379_ap 


the nucleus was also proposed by [Basilevsky & Kell^ ( |2007| ) 
for comets 19P/Borrelly, 81P, and 9P. 

The tensile, shear, and compressive strengths of the cometary 
material have been estimated by several methods, including the 
Deep Impact experiment, comet breakup observations, labora¬ 
tory experiments, and theoretical modeling (Table l). |Biele et al.| 
( [2009| ) compiled and discussed these different estimates in an 
excellent review paper. From the Deep Impact experiment, the 
shear strength was estimated to be <65 Pa by | A’Hearn et al.| 
( [2005 1), but might be any value b etween 0 and 12 kPa according 
to [Holsapple & Housenj ( [2007 j ); this is an uncertainty of three 
orders of magnitude. From comet breakup, the tensil e strength 
was estimated to be 5 Pa by [Asphaug & Benzj ( |1996| ) from the 
encounter of comet Shoemaker-Levy 9 with Jupiter, to 100 Pa 
for a sun-grazing comet with a radius of 1 k m ([Klinger et al" 


[1989[ ). [Toth & Lisse[ ( [2006[ ) and [Davidsson[ ( [2001[ ) estimated 
that a tensile strength of <100 Pa and 1-53 Pa, respectively, is 
sufficient to keep a comet nucleus stable against its rotational 
breakup. From laboratory experiments, the compressive strength 
of cometary material analogs (water ice and dust mixture) was 
estimated to be between 20kPa and 1 MPa ( [Jessberger & Kot-[ 
[thaus| 198^[Bar-Nun et al.|20'07] ), while the tensile strength was 
estimated to be 200-1100 Pa for homogeneous Si02 dust sam¬ 
ples ( Blum et al.[2006 ) and down to 1 Pa for dust-aggregate sam¬ 
ples ( Blum et al.|2014[ ). From modeling, the tensile strength of 
fluffy silicate dust/ice material was estimated to be 270 Pa by 
[Greenberg et al.[ ( p^95] ) and to be 5000Pa by [Biele et al.[ ( [2009] ) 
for the same material made of water ice alone. Again from mod¬ 
eling, the compressive strength was estimated to be 6500 Pa for 
porous icy grains ( [Sirono & Greenberg |2000[ ). 

Depending on whether the material is consolidated or un¬ 
consolidated, these strength estimates (Table 1) depend, or not, 
on the scale at which they were measured. While the strength 
of unconsolidated material might be scale invariant, that of con¬ 
solidated material follows a typical d~^ power law, where d is 
the scale and q is the power exponent, with q~0.6 for water ice 
( [Petrovic|2003[ ). From the kilometer to the millimeter scale, the 
strength of consolidated material can thus change by three orders 
of magnitude, still less than the above ranges, which cover up to 
five orders of magnitude. The different estimates clearly are a 
priori difficult to reconcile with each other, and large uncertain¬ 
ties remain on the tensile, shear, and compressive strengths of 
the cometary material. 

Section 2 presents the data, shape models, and methods used 
in this paper. Section 3 discusses the link between gravitational 
slopes and surface morphologies for different types of terrains 
and regions on the nucleus. In Sect. 4 we estimate the tensile, 
shear, and compressive strengths of the cometary material. Dis¬ 
cussions and conclusions are presented in Sect. 5. 


2. Data, shape models, and gravitational slopes 

All the images shown in this paper were acquired with the Nar¬ 
row Angle Camera (NAG) of the Optical, Spectroscopic and 
Infrared Remote I maging System (OSIRIS) onboard Rosetta 
( Keller et~aL][2007[ ) since August 2014. Their spatial resolution 
varies between 18cmpix“^ and 1.8mpix“^ 

For this work we used two different shape models of the nu¬ 
cleus of 67P, computed from OSIRIS images. The first shape 
model was computed by [Jorda et ak] p015[ ) using the stereo- 
photoclinometry te chnique (SPC). The se cond shape model 
was computed by [Preusker et al.[ ( [2Q15[ ) using the stereo- 
photogrammetry technique (SPG). The SPC shape model was 
resampled to match the resolution of the SPG shape model. The 
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Digital Terrain Models (DTMs) extracted from the SPC and SPG 
shape models have a horizontal sampling of 2 m and a typical 
vertical accuracy at the decimeter scale. 

We computed the local gravitational slopes for these two 
shape models, including the effects of the nucleus rotation and 
assuming a uniform density inside the nucleus. Details on the 
method are provided in [ Jorda et al.[p012) . The error on the grav¬ 
itational slope is estimated to be 5°. When not specified, the term 
slope in this paper always refers to the gravitational slope. 


3. Relation between gravitational slopes and 
surface morphologies 


The geomorphology of the nucleus surface is diverse and con¬ 
strained by several processes related to gravity and cometary ac¬ 
tivity ( [Sierks et al.|2015[ [Thomas et al.|2015[ ). The link between 
the different types of terrains (smooth, hummocky, consolidated 
material, dust covered, with or without boulders, etc.) and their 
gravitational slope is important to better constrain the processes 
in play and the nature of the cometary material. In this section 
we study this relationship for five regions on the nucleus that are 
representative o f the different morpholog ies observed on the sur¬ 
face. We refer to [El-Maarry et al.[ ( |2015a ) for the definition of the 
regions. 


1. The Imhotep region (Fig.[^ - This region presents a wide 
variety of terrains and morphologies. The most remarkable 
ones are the smooth terrains, the largest of which extend over 
0.8 km^ , and the roundish features observed near the gravi¬ 
tational low of the region. 

2. The Ash region (Fig.[^ - This region is mostly covered by 
dust that is spatially unresolved at the decimeter scale. It 
shows several debris fields that are made of boulders, which 
are located at the feet of steep walls exposing a consolidated 
material. A large depression of 370 m width is visible on the 
left side (noted A in Fig.[^. 

3. The Seth region (Fig.[^ - This region is dominated by cir¬ 
cular depressions, most of them being accumulation basins 
with an opening toward a lower basin. Debris accumulates at 
the feet of the steep walls of basins. A large (-200 m) and 
deep pit dominates the bottom right part of the region (noted 
A in Fig.[^. 

4. The Hathor region (Fig.|^ - This region is dominated by 
cliffs with a maximum height of -900 m, exposing a consol¬ 
idated material. Many anisotropies (e.g., fractures) are visi¬ 
ble, organized in two main directions roughly perpendicular 
( [Thomas et al.|2015[ ). 

5. The Agilkia region (Fig.|^ - Agilkia is a large “super region” 
corresponding to the nominal landing site and includes the 
regions of Hatmehit, Ma’at, Nut, and Maftet. The Agilkia re¬ 
gion is dominated by a very large depression of 800 m width 
on its right side (noted A in Fig.[^. This region shows a vari¬ 
ety of exposed consolidated material, dust-covered material, 
and boulder fields. 


The histogram of gravitational slope angles computed for 
each region is shown in Fig.|^ Each region has a unique dis¬ 
tribution of gravitational slopes, related to its unique topogra¬ 
phy. There are similarities between the different distributions, 
however. All regions cover a wide range of slopes from 0° to 
almost 100°. All distributions show a peak (Eig.[^ right panel). 
The peak of the distribution varies between 2° for Imhotep and 
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64° for Hathor. The peak value is a good indication of the overall 
flatness of the terrain, lower values indicating a flatter terrain on 
average. Secondary peaks, bumps, or shoulders are also visible; 
they are related to steep walls (e.g., bump around 55° on Seth 
and Imhotep) or smooth terrains (e.g., shoulder around 20° on 
Imhotep and Agilkia). Slopes in excess of 90° are indicative of 
overhangs and are only visible on the Hathor cliffs; overhangs 
cover 0.1% of the total Hathor area. 

Figurej^ allows a direct comparison between the SPC and 
SPG gravitational slopes. This is important to better estimate the 
error on the slope. The overall shape of the SPC and SPG slopes 
distributions is similar. In particular, the peaks are at the same 
position. The notable exception is Hathor, with a peak at 64° 
for the SPC slopes and 70° for the SPG slopes. This leads to a 
more general comment, which is that the slope distributions of 
SPC and SPG differ for slopes steeper than 60°. The SPC shape 
model tends to systematically underestimate the fraction of the 
surface with steep slopes compared to the SPG shape model. 
This exercise shows that slopes are robust up to 60°, with an 
error of +5° as given in Sect. 2, but larger uncertainties affect 
slopes above 60°, with an error up to 20° in some cases. 

Figures[^-[^illustrate that all regions show a similar general 
pattern in terms of the relation between gravitational slope and 
terrain morphology. A sketch of this general pattern is shown in 
Fig-i 

- Low-slope terrains. Terrains with slopes in the range 0 - 20° 
are covered by a spatially unresolved material, that is, a ma¬ 
terial made of particles smaller than 20 cm that we call the 
fine material in this paper. A few large isolated boulders are 
visible on these terrains, with a typical size larger than 10 m. 


“ Intermediate-slope terrains. Terrains with slopes in the range 
20-45° are mainly fallen consolidated materials and debris 
fields, with numerous intermediate-size boulders from <1 m 
to 10 m for the majority of them (see also |Pajola et al.|20T5] ). 
These terrains are covered by a dust deposit of variable thick¬ 
ness, which partially hides some boulders. Most of these ter¬ 
rains are located at the feet of high-slope terrains. 

- High-slope terrains. Terrains with slopes in the range 45 - 
90° are cliffs that expose a consolidated material and do not 
show boulders or a fine material. These terrains probably 
show the bare nucleus. 

This general pattern is very similar to what we And on Earth, 
particularly in young mountains like the Alps, where boulder 
fields are frequently observed at the feet of cliffs; theses fields 
are a result of cliff collapse. The sublimation of ices triggers ero¬ 
sion, most likely exacerbated by fractures; gravity controls the 
collapse like it does on Earth. 

The presence of large boulders on slopes lower than 20° is 
intriguing (yellow circles in Figs. 1 These boulders are usu¬ 
ally isolated, far from high-slope terrains, and are large, tens 
of meters. They are too large to be lifted by gas drag resulting 
from regular cometary activity, this m echanism only applies to 
boulders smaller than the meter scale ([Groussin & Lamy|2003[ 
lA’Heam et aq|MTT| [Kdky et al.||20T^Gundlach et al.||20I5| ). 

They could instead be air falls from outburst events, which are 
limited in time, but much stronger in terms of released energy. 
Such events are sporadic and spatially localized, however; they 
may not be frequent enough to explain the presence of large 
boulders in all the regions we studied. An alternative solution 
is that large boulders are leftovers from previous basins and de¬ 
pression edges, when they were smaller and less eroded. 


The transition between intermediate- and high-slope terrains 
is sharp. This indicates that high-slope terrains have a steeper 
slope than the angle of repose, which can be es timated to b e 
45+5°, a typical value for gravel on Earth (|Julien|I995 ). 




repose 


The transition between low- and intermediate-slope terrains is 
softer, with several examples of boulder fields on slopes lower 
than 20° (white circles in Fig.[^ and of smooth terrains on slopes 
steeper than 20° (red circles in Figs. [I]-0. 

Boulder fields that end on slopes lower than 20° on Seth 
(white circles in Fig.[^ could result from a progressive erod¬ 
ing and degradation process of the nucleus. These boulder fields 
were probably rock falls, at the feet of the previous emplacement 
of cliffs, as the above isolated boulders. Following a slow degra¬ 
dation process that involves fractal fragmentation into small 
pieces and progressive covering by dust deposits, t hey now ap- 
pear on low-slope terrains and are partially hidden ( [Pajola et"ar 


|20I5| ). The same degradation process can explain that a fine 
material of variable thickness partially hides boulder fields on 
intermediate-slope terrains. 

Smooth terrains on slopes steeper than 20° (red circles in 
Figs.[T]-|g could be ancient boulder fields that were able to retain 
the dust deposit on a slope steeper than that of typical smooth 
terrains. The origin of the dust could be deposits from regular 
cometary activity or products of the degradation process of cliffs 
and boulders. 


4. Tensile, shear, and compressive strengths of the 
cometary material 

4.1. Definition of strengths 

There are three types of strength for a given material: the tensile 
strength ctj, the shear strength , and the compressive strength 
(Tc- They define the ability of a material to withstand mechanical 
constraints. Usually, ctt < crs < crc. Depending on the nature 
of the material, consolidated or unconsolidated, the strength de¬ 
pends, or not, on the scale at which it is measured (Sect. 1). A 
review on sca ling effects on structural strength can be found in 
|Bazant| ( [l999| ). 

As explained in the introduction, large uncertainties remain 
in the values of ctt, crs , and crc for the come tary material, some - 
times of several orders of magnitude (e.g., |Biele et ^|2009| ). 
We here intend to provide additional constraints on the strengths 
of the cometary material; more precisely, the tensile strength of 
overhangs and collapsed structures, the shear strength of fine ma¬ 
terials, boulders and Hathor cliffs, and the compressive strength 
of fine materials and consolidated materials. 


4.2. Tensile strength of overhangs and collapsed structures 


The tensile strength can be estimated from overhangs. From sim¬ 
ple mechanics (e.g., |Tokashiki & Aydan|2010| ), the failure of an 
overhang of rectangular shape due to bending will occur if the 
following condition is fulfilled: 


ctt < By — 
^ H 


( 1 ) 


where y (N m“^) is the unit weight, L (m) the length of the over¬ 
hang, and H (m) its height as defined in Fig.[^ Equation ([^, 
which is independent of the width X of the overhang, can be 
translated into 


(Tt < 3pg 


H 


( 2 ) 
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where p = 470 kg is the density of the material (Sierks et al. 
2015) and g is the local gravity. 

Several examples of overhangs or collapsed structures are 
shown in Figs.p^andpT] Figurep^ shows two large structures 
that collapsed, indicating that the tensile strength was exceeded. 
Using the DTM of the area, we estimated the length and height of 
the largest structure to be L = 100 m and = 30 m. From Eq. ^ 
with g = 2x 10“^ m s“^ , this gives an upper limit for the tensile 
strength ctj < 94 Pa. In Fig.[TT] several overhangs are visible 
with an estimated length of 10 m and a height of 5 m. The pres¬ 
ence of mass wasting in the form of boulders at the feet of these 
overhangs indicates that they are close to breaking and thus good 
estimates of the tensile strength. This argument is also supported 
by the fact that some of these boulders have a size (~10m) sim¬ 
ilar to that of the overhangs themselves. From Eq. Q we derive 
(Tt = 5.6 Pa for these overhangs. 

This determination of the tensile strength is not accurate, 
mainly because of the uncertainty in the geometry (length and 
height). Taking these uncertainties into account, the tensile 
strength of small overhangs (~10m) is most likely in the range 
3-15 Pa, and the tensile strength of large collapsed structures 
(~100m) is lower than 150 Pa. 


4.3. Shear strength of fine materials, boulders, and Hathor 
cliffs 

The shear strength can be calculated with Eq. Q, where 6 is the 
slope angle on which the boulder is located, m (kg) is the mass 
of the considered boulder, and A (m^) is the contact area of the 
boulder with the terrain underneath. The shear strength results 
from friction and cohesion, and we cannot separate these two 
physical quantities in this study. 


crs 


mg sin 6 


(3) 


For a boulder of radius r, Eq. ([^ translates into Eq. Q with A = 
7r(rcos(p)^ (Fig. 12). 


crs = 


^Trr^pg sin 0 4rpgsin0 


7i(r cos ipY 


3 cos^ ip 


(4) 


The shear strength is constrained by the largest boulders on 
the highest gravitational slopes, more precisely, those with the 
highes t rsin^ value. From Fig.[^ adapted from | Auger et al.| 
( |2015| ), the highest value for rsin^ is 5.2 m and corresponds to 
a boulder with radius r = 11.5m located on a gravitational slope 
6 = 26.8°. Assuming an area of contact of 1 % of the total boul¬ 
der surface, defined with p = 80° in Eq. Q, we obtain a shear 
strength crs = 22 Pa. Because of the uncertainty on the radius 
(2 m), slope angle (5°), and the fact that the area of contact may 
be larger than I % (p < 80°), the shear strength of the boulders 
and fine material on which they stand is most likely in the range 
4 - 30 Pa. Figure[^also illustrates that there are indeed no boul¬ 
der on slopes steeper than the angle of repose ^repose =45 + 5° 
(Sect. 3). 

A lower limit of the shear strength of the Hathor cliffs is pro¬ 
vided by the lateral pressure at the bottom of these cliffs, given 
by Eq. 0, where h is the height of the cliff. 

os > pgh(l - sin6»repose) (5) 


For a maximum height h = 900 m of the Hathor cliffs, this gives 
CTs > 30 Pa. This value is an approximation since the cliff does 
not have a slope of 90° from top to bottom, meaning that it is not 
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perfectly vertical, and the gravity changes from top to bottom. 
Nevertheless, it shows that due to the low gravity, a low strength 
is sufficient to withstand a high cliff. 


4.4. Compressive strength of fine and consolidated materials 

The compressive strength can be estimated from the footprints 
left by the lander Philae after it bounced on the nucleus surface 
at the first nominal landing site (Fig. [g. We can derive the com¬ 
pressive strength with Eq. Q, where F (N) is the forced applied 
to a surface of area A (m^), m (kg) is the mass of Philae, v (m s“^) 
is the lander impact velocity, and d (m) is the depth at which the 
lander penetrated the surface material. 



1 2 
dA 


( 6 ) 


The depth of the largest footprint, which has a diameter of 
~5 pix or 150 cm, was estimated to be 20 cm using a shape from 
shading analysis. This footprint was created by Philae when the 
first leg touched the ground; we assume here that the three legs 
did not touch the ground simultaneously. With a mass of 100 kg 
for the lander, an impact velocity v = 1 m s“^ and a contact area 
of 0.016 m^ per foot (each foot is made of two disks of 10 cm 
diameter), we derive a compressive strength of 15.6 kPa for the 
surface material. This value is an upper limit since the presence 
of footprints indicates that the impact process exceeded the com¬ 
pressive strength of the surface material. Moreover, as we dis¬ 
cuss in Sect. 5.4, the penetration of Philae of the surface mate¬ 
rials was blocked by a hard layer under the surface, which also 
leads to overestimating the compressive strength. 

A more reliable estimate can be given by taking into ac¬ 
count mechanical considerations. In practice, the compressive 
strength of a consolidated material is always larger than its ten¬ 
sile strength by typically one order of magnitude (e.g., rocks like 
sandstone), so that we can reasonably assume crc ~ IOctj^ ( |Sheo- 
rey|1997| ). This only applies to consolidated materials, however, 
such as overhangs or collapsed structures, but not to a uncon¬ 
solidated material like the fine material on the surface ( |Blum &| 
|Schrapler|2004t[Blum et al.|2006| ). Prom the tensile strength cal¬ 
culated in Sect. 4.2, we therefore derive a compressive strength 
of the consolidated material that most likely is in the range 30 - 
150 Pa, with an upper limit of 1.5 kPa; this is well below the 
upper limit of 15.6 kPa derived for the fine surface materials. 


5. Discussions and conciusions 

5.1. Summary and comparison with other measurements 


Values for the tensile, shear, and compressive strengths are 
summarized in Table[^ Overall, the strengths are low and the 
cometary material can be considered as weak. Once scaled to 
the meter scale for comparison, our results agree well with most 
estimates from other authors who used observations, laboratory 
experiments, or modeling, in particular for the tensile strength. 
This scaling should be taken with caution, however, since con¬ 
solidated and unconsolidated materials may not follow the same 
scaling law. Our estimate for the shear s trength is at the lower 
end of the possible range mentioned by |Holsapple & Hous^ 
( |2007|). Concerning the com pressive strength, the discrepancy 
with |Jessberger & Kotthaus|(l989| ) results from the fact that they 
measured the compressive strength of a mixture made of liquid 
water and dust with a high fraction of water (>80%). This mix¬ 
ture is cooled from room temperature to low temperature (253 K, 




























Groussin, Jorda, Auger et al.: Gravitational slopes, geomorphology and material strengths of the nucleus of comet 67P 


233 K, or 123 K, depending on their test case), and water ice 
forms as compact hexagonal ice, known to be hard (crc = 5- 
25 MPa; |Petrovic|2003) . 

5.2. Strength-to-gravity ratio 

The strengths are low on 67P, but the gravity is low as well. The 
ratio between the tensile strength at the 1 m scale (7 - 34 Pa) and 
gravity (2 x m s“^) is 35 -170 x 10^ Pa s^ m“^ on 67R On 
Earth, weak rocks like siltstone, with a density of ~2600 kg m“^ 
and a porosity of 21-41 %, have a typical tensile strength of 
0.5 MPa. For this type of weak rocks, the ratio between the 
tensile strength and Earth’s gravity is 51 x 10^ Pas^ m“^. This 
remarkable agreement between the strength-to-gravity ratio on 
67P and on Earth explains why the general pattern between grav¬ 
itational slopes and surface morphologies looks “familiar” to 
our eyes, from high-slope terrains (cliffs) to intermediate-slope 
terrains (mass wasting as boulder fields) and low-slope terrains 
(fine material). To first order, gravity shapes terrains in a similar 
way on 67P and on Earth. 


5.3. Low strengths and surface morphologies related to 
activity 


The low strength of the cometary material may help to under¬ 
stand some of the surface morphologies, in particular those im¬ 
plying the rising up of a gas bubble, such as the roundi sh features 
of 67P (Imhotep region; [Auger et al.|2015| ) and 19P ( [Brownie^ 
jet al.|2Q0^ , or the depressions adjacent to smooth terrains on the 
nucleus of 9P ( [Belton & Melosh||2009| ). The saturated vapour 
pressure of the gas inside the nucleus exceeds 1 kPa when the 


temperature exceeds 55 K for CO and 155 K for CO 2 (Prialnik 


[et al.|2004| ). This means that the CO and CO 2 gas pressure may 
locally exceed the compressive strength of the cometary material 
and the gas may push the surrounding materials while it rises to 
the surface, creating channels inside the nucleus. This process re¬ 
mains speculative, however, since it has a severe limitation: how 
would solar energy penetrate meters below the su rface, despite 


the low thermal inertia of 10-50 JK ^ m ^ s ([Gulkis et al. 


2015)? We currently do not have the answer to this critical ques¬ 


tion, but the presence of fractures on the surface ( [Thomas et al. 
2015[[ErMaarry et al. 12015b] ) may help to solve this issue. 


5.4. Philae observations 

A key question is how to reconcile our low strength estimates 
with the observation made by the Philae MUPUS experiment, 
which suggests much larger strengths of several MPa since MU¬ 
PUS was not able to penetrate the surface material below a few 
centimeters ( [Spohn et al.|2014[ ). The solution resides in the hard 
la yer of water ice under t he surface dust deposit. As explained 
in [Pommerol et al. ( 2015[) , laboratory experiments with analogs 
such as KQSI ( [Grun|1991[ ) have shown that a hard layer of water 
ice can be produced by sublimation/redeposition cycles and/or 
sintering of water ice close to the surface. This layer has an es¬ 
timated thickness of a few centimeters ( Pommerol et al.|[2015| ) 
to several meters ( [Kossacki et al.|[MT5 ) and has been heavily 
processed, so that it is now hard with strengths up to 1 MPa 
as measur ed by [Jessberger & Kotthaus| ( [1989[ ), [Seiferlin et al.[ 
( [1995[ ), or |Kochan et al.| ( [1989 ). This layer, which is not repre¬ 
sentative of the bulk nucleus material, provides a simple expla¬ 
nation for the MUPUS observations. The presence of this layer 
is compatible with our results on strength since we derived i) the 


tensile strength from large overhangs and collapsed structures, 
for which this layer is likely negligible, otherwise the measured 
tensile strength would be much larger, ii) the shear strength from 
boulders on the surface, for which this buried layer should not 
matter, and iii) the compressive strength from the tensile strength 
using mechanical considerations, for which the same justifica¬ 
tion can be used. 


5.5. What is representative of the bulk nucleus material? 


We derived strengths for different materials in different locations 
on the nucleus, and the question naturally arises whether these 
materials are representative of the bulk nucleus material. While 
it is not possible to give a firm positive answer to this question, 
we provide here some arguments supporting this idea. 

First, it is very satisfying to see in Table 1 that strength esti¬ 
mates derived for different materials at different scales and with 
different methods are mostly consistent with each other, and this 
from unconsolidated material with modeling at the microme¬ 
ter scale to consolidated material with observations at the meter 
scale. When translated to the same scale, dust aggregates, peb¬ 
bles, boulders, or larger structures have similar strengths. Since 
they are the building blocks of the nucleus, their mechanical 
properties should be representative of that of the bulk cometary 
material. 


Second, we determined the tensile strength for some large 
collapsed structures, with a thickness of up to 30 m. These large 
structures should contain a significant fraction of material that 
has not been, or has only partially been, affected by activity pro¬ 
cesses related to insolatio n As a result of the l ow thermal inertia 
of 10-50 JK“^ m“^ s“^-^ ( Gulkis et al.[[2015 1, the thermal heat 
wave only penetrates the nucleus by a few meters (at most) (e.g., 
Groussin et al.[[^13| ). This is confirmed by modeling, which 


shows that only the first 5 meters are indeed affected by inso¬ 
lation ( [Prialnik et al. [2004} . Complex models also show that wa¬ 
ter ice is buried under a thin la yer (cm) of material depleted in 
volatiles ( [Skorov & Blum[20T^ . The presence of water ice close 
to the surface (~ 1 m) was also observed by the Deep Impact ex¬ 
periment ( [Sunshine et al.|2007] ). All these arguments support the 
idea that activity processes only affect the first meters of the nu¬ 
cleus. The large 30 m thick structures for which we derived the 
tensile strength should therefore be a good proxy for the bulk 
material of the nucleus, even if they are not fully representative 
of it. We furthermore add that it is currently unclear whether the 
alteration of the first meters would make the material weaker 
through forming fractures, or stronger through forming a hard 
sintered layer. 


Third, it i s important to add tha t, as demonstrated by [Toth & 
|Lisse| ( [2006] ), [Toth & Lisse[ ( |2010| ), and ( [Davidsson|2001[ ), even 
a low tensile strength of a few hundred Pa is amply sufficient 
to prevent the rotational breakup of the nucleus of 67P with its 
current rotation period of 12.4 h. 


Finally, although we currently have no evidence for it, we 
cannot exclude heterogeneities in strengths over the nucleus. For 
example, the Hathor region, for which we determined a shear 
strength >1770 Pa at the 1 m scale, may be stronger than some 
small overhangs. Addressing this problem requires more exten¬ 
sive investigations and is beyond the scope of this paper. For 
now, we can only argue that the general pattern of Fig.[^ which 
is directly linked to the material strengths, is valid across the en¬ 
tire nucleus. 
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5.6. Constraining the origin of 67P 

It is interesting to compare the compressive strength to the 
pressure resulting from gravity inside the nucleus. Assuming a 
constant density, the pressure P inside the nucleus is given by 
Eq.Q, where R is the nucleus radius and r is the distance from 
the center. 


P{r) = -nGp^R^ 


1 




(7) 


Figure[^ shows the pressure inside the nucleus as a function of 
depth for different initial radii of the nucleus after its complete 
accretion and before it entered the inner solar system. Depend¬ 
ing on past evolution scenarios, the nucleus of 67P was initially 
larger by a few hundred meters to a few kilometers ( [Groussin] 
|et al.||2007] ). The nucleus has shrunk as a result of erosion and 
now has a radius of ~2 km (red line in Fig. [T^. For an initial 
radius of 2.2 - 3.0 km, slightly larger than today, the pressure in¬ 
side the nucleus exceeded the compressive strength of the layers 
we see today (30-150 Pa). This means that if the nucleus was 
formed by the gentle accretion of pristine materials (see below), 
their subsequent compaction by gravity inside the nucleus was 
sufficient to explain the compressive strength we measure today. 

Continuing with this idea, one may ask whether compaction 
by gravity inside the nucleus was also sufficient to form the lay¬ 
ers we see today on the nucleus ( [Massironi et al.|[2015| ). On 
Earth, the diagenesis of rocks starts when they are stressed to 
more than ten times their compressive strength. If the same 
mechanism applies to 67P, the primordial material that accreted 
to form the nucleus had a compressive strength at least ten times 
lower than the one we measure today to initiate diagenesis, that 
is, as low as 3- 15 Pa. In these conditions, compression could 
have contributed to the for mation of layers. _ 

Following the model of |Skorov & Blum| ( |2Q12| ), the tensile 
strength of dust aggregates formed by gravitational instability 
typically is 1 Pa. This value is remarkably consistent with our re¬ 
sults (Table 1), particularly if we take into account a compaction 
scenario where the primordial material was even weaker. Our 
results therefore tend to favor a formation of comets by pebble 
accretion in a region of higher concentration of particles such as 
vortices ([Johansen et al.|20f4[ [Barge & Sommeria||1995[|Blum| 


et al.|2014| ), which implies a gentle formation process by accre 

tion at low velocity, on the order of 1 m s“^ or lower. In con 
trast, the hierarchical accretion model ( |Weidenschilling|[2004| ) 
with velocities up to 50ms“^ for particles larger than 1 m, and 
the collisional scenario between two large bodies of tens of km 
or more ( [Davis & Farinella|1997| ) with an internal compression 
by gravity larger than lOkPa, although not excluded, are less fa¬ 
vored. 

Finally, the low tensile strength indicates that the nucleus of 
67P very likely did not experience a global melting/freezing of 
its water content, neither during the accretion stage by radio¬ 
genic heating ( jPodolak & Prialnik|2000| ), nor more recently by 
the exothermic amorphous-to-crystalline water ice reaction ( jPri-j 
alnik et ar]|2004| ). Indeed, if a global melting/freezing had oc 


curred, the nucleus would be much harder, from 10 kPa ( [Miles &| 
Faillace|2012[ ) to 1 MPa (e.g., [Seiferlin et al.|1995[), and denser 


with a bulk density higher than 1000 kg m ^ ([Miles & Faillace 

[MT2| ). 
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Table 1. Summary of the tensile strength (cr^), shear strength ( 07 - 5 ), and compressive strength (ctc) for cometary materials, including our own 
estimates. The list is not exhaustive, but representative of values found in the literature. We scaled each value to the meter scale using a power law 
with an exponent equal to -0.6 (Sect. 1). This scaling should be taken with caution since consolidated and unconsolidated materials may not follow 
the same scaling law. 


Reference 



Method 

Value (Pa) 

Scale 

Value (Pa) 
at 1 m scale 

Tensile stren 

Klinger et al. 

gth (cTj.) 

(19891 

Observations (sungrazing comets) 

100 

Im 

100 

Davidsson (|2 

DTJT 

1 

Observations (rotational breakup) 

1-53 

Im 

1-53 

Asphaug <K Benz 

(19961 

Observations (D/1993 F2 SL9) 

5 

Im 

5 

Blum et al. (|2006 


Laboratory experiments 

200-1 100 

1 cm 

13-69 

Bar-Nun et al. (|2t 

)07l 

Laboratory experiments 

2000-4000 

lOOyum 

8-16 

Blum et al.ldl^ 

[DPT 


Laboratory experiments 

1 

1 mm 

0.02 

Biele et al.|(|2 

009| 

1 

1995ji 

Modelling 

5 000 

lOyum 

5 

Cireenberg et 

al. 1 

Modelling 

270 

lOyum 

0.3 

This work 



Observations (67P, Rosetta) 

- Collapsed structures 

- Overhangs 

<150 

3-15 

30 m 

5 m 

<1150 
8-39 


Shear strength (cr^) 

Holsapple & Housen (2007) 
A’Hearn et al. (|20(b|i - 

This work 

Observations (9P, Deep Impact) 
Observations (9P, Deep Impact) 

Observations (67P, Rosetta) 

- Hathor cliffs 

- Fine surface materials and boulders 

0-12000 

65 

>30 

4-30 

Im 

Im 

900 m 
Im 

0-12000 

65 

>1770 

4-30 

Compressive strength (ctc) 
Jessberger & Kotthaus (19891 

Laboratory experiments 

30000-1000000 

1 cm 

1890-63100 

Bar-Nun et al. (1200 / 1) 

Guttler et al. (|2D09|l 

Laboratory experiments 

20000 

lOOyum 

80 

Laboratory experiments 

<1000 

1 mm 

16 

Blum et al. (|z0l4|) 

Laboratory experiments 

15 

1 mm 

0.2 

Sirono & ureenberg (2000) 

Modelling 

6500 

lOyum 

7 

This work 

Observations (67P, Rosetta) 

- Fine surface materials 

<15 600 

10 cm 

<3 920 


- Collapsed structures 

<1500 

30 m 

<11600 


- Overhangs 

30-150 

5m 

79-394 
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Groussin, Jorda, Auger et al.: Gravitational slopes, geomorphology and material strengths of the nucleus of comet 67P 



Fig. 1. Left panel: view of the Imhotep region. The letter A indicates the region of roundish features, observed near the gravitational low of the 
region. Image NAC_2014-08-25T23.12.54 (spatial resolution: 0.94 mpix"^). Right panel: gravitational slopes of the Imhotep region, computed 
from the SPC shape model ( [Jorda et al.|2015| ), superimposed on the background image of the left panel. Blue corresponds to terrains with slope 
angles between 0° and 20°, yellow to terrains with slope angles between 20° and 45° , and red to terrains with slope angles between 45° and 90°. 
Yellow circles show examples of large and isolated boulders on slopes lower than 15°. White circles show examples of boulder fields that end on 
slopes lower than 15°. 



Fig. 2. Left panel: view of the Ash region. The letter A indicates the largest depression of the region. Image NA C_2014-08-07Tl 8.20.34 (spatial 
resolution: 1.5 mpix"^). Right panel: gravitational slopes of the Ash region, computed from the SPC shape model ( jJorda et ah 2015 1, superimposed 
on the background image of the left panel. Blue corresponds to terrains with slope angles between 0° and 20°, yellow to terrains with slope angles 
between 20° and 45° , and red to terrains with slope angles between 45° and 90°. Yellow circles show examples of large and isolated boulders on 
slopes lower than 15°. Red circles show examples of smooth terrains on slopes in the range 25 - 30°. 
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Fig. 3. Left panel: view of the Seth region. The letter A indicates the largest pit of the region. Image NAC_ 2014-08-16T10.59 .16 (spatial resolution: 
l.Smpix"^). Right panel: gravitational slopes of the Seth region, computed from the SPC shape model ( Jorda et al.|2015| , superimposed on the 
background image of the left panel. Blue corresponds to terrains with slope angles between 0° and 20°, yellow to terrains with slope angles between 
20° and 45° , and red to terrains with slope angles between 45° and 90°. Yellow circles show examples of large and isolated boulders on slopes 
lower than 15°. Red circles show examples of smooth terrains on slopes in the range 25-30°. White circles show examples of boulder fields that 
end on slopes lower than 15°. 
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Fig. 4. Upper panel: View of the Hathor region. Image NAC_2014-08-07T20.20.34 (spatial resolution: l.Smpix"^). The red circle shows an 
example of a smooth terrain on slopes in the range 25 - 30°. Green circles show examples of overhangs with slopes >90°. Lower panel: Ext racted 
Digital Ter rain Model (left) of the Hathor region with the corresponding gravitational slopes (right) computed from the SPC shape model ( |Jorda| 
|et al.|2015] ). For technical reasons related to the DTM extraction, we were not able to produce a figure for Hathor like Figs.[^[^[^ or|^ 
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Fig. 5. Left panel: view of the Agilkia region. The letter A indicates the largest depression of the region. Image NAC_201 4-08-25T15.42.54 
(spatial resolution: 0.94 mpix"^). Right panel: gravitational slopes of the Agilkia region, computed from the SPC shape model ( |Jorda et ^ 20151, 
superimposed on the background image of the left panel. Blue corresponds to terrains with slope angles between 0° and 20°, yellow to terrains 
with slope angles between 20° and 45° , and red to terrains with slope angles between 45° and 90°. Yellow circles show examples of large and 
isolated boulders on slopes lower than 15°. The red circle shows an example of a smooth terrain on slopes in the range 25 - 30°. 




Gravitational slope angle (deg) 


Fig. 6. Histogram of gravitational slope angles for the five regions, with one color per region, in logarithmic scale (left) and linear scale {right) for 
the Y axis. The values are normalized so that the sum of all Y values equals 1. Each region has a unique distribution of gravitational slopes, related 
to its unique topography. All re gions cover a widee range of slopes from 0° to almost 100°. The gravitational slopes were derived from the SPC 
shape model ( |Jorda et al.|2015| ). 
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0 20 40 60 80 100 

Gravitational slope angle (deg) 



Gravitational slope angle (deg) 


Fig. 7. Histograms of gravitational slope angles for the five regions, derived from the SPC shape model (solid line) and from the SPG shape model 
(dashed line). The black arrow on Hathor indicates a bump on the SPG slopes distribution that must be ignored since it results from the fact that the 
SPG DTM includes a large part of the Hapi region that is absent from the SPC DTM. The Hapi region, which is smooth and fiat, thus introduces a 
bias in the distribution for low-slope angles. 
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Fig. 8. Sketch of the general pattern in terms of relation between gravitational slopes and terrain morphologies: i) low-slope terrains (0-20°) 
are covered by fine material and contain a few large (>10m) and isolated boulders, ii) intermediate-slope terrains (20-45°) are mainly fallen 
consolidated materials and debris fields, with numerous intermediate-size boulders from <1 m to 10 m for the majority of them, and iii) high- 
slope terrains (45-90°) are cliffs that expose a consolidated material and do not show boulders or fine material. The borders between low- and 
intermediate-slope terrains are not always sharp (see text for details). 



Fig. 9. Scheme of an overhang with rectangular shape L x X x H. 
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Fig. 10. Examples of large structures that already collapsed, with a length L of 100 m or 50 m and a height // of 30 m (determined from the Digital 
Terrain Model). These structures are located at the northern border of the Imhotep region. Image NAC_20I4-09-I6TI5.44.07 (spatial resolution: 
0.52 mpix"^). 


Article number, page 15 of 18 







ASlA proofs: manuscript no. 26379_ap 



Fig. 11. Red/blue anaglyph of overhangs at a high spatial resolution of IScmpix"^ in the Maftet region. Several overhangs are visible in this 
anaglyph; they are indicated by white arrows, all with a length of ~10m and a height of ~5 m. Image NAC_2014-10-19T13.18.55 (spatial resolu¬ 
tion: 0.18 mpix"^). 
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Fig. 12. Scheme of a boulder lying on the surface. The radius of the boulder is r, and its contact area with the surface is defined by the angle (p. 
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Fig. 13, Lef t panel: size (radius r) of 2207 boulders of the Imhotep region as a function of their gravitational slope angle 6 (adapted from [Auger 
jet al.|2015|). The black horizontal lines correspond to the angle of repose ^repose = 45 + 5°; there are no boulder on slopes larger than the angle of 
repose. Right panel: value r sin 6 for each boulder, as a function of its radius. Both panels: for robustness, we rejected the upper 1% of the r sin 6 
distribution (blue points). For the remaining 99%, the boulder with the highest value r sin 6 is highlighted in green and corresponds to r = 11.5 m 
and 6 - 26.8°, as indicated by the green horizontal and vertical dashed lines. 
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Before After After (zoom) 
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Fig. 14. Images of the first nominal landing site, before and after the Philae touchdown. The footprint of the three landing feet are clearly visible. 
Image NAC_2014-11-12T15.18.52 and NAC_2014-11-12T15.43.51 (spatial resolution: 0.31 mpix-'). 
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Fig. 15. Pressure versus depth inside the nucleus, for different initial sizes of the nucleus (radius R-2.2, 3.0, 5.0 or 7.0 km). The pressure increases 
with depth. The blue range corresponds to the best value of 30-150 Pa for the compressive strength of the surface materials we see today (30- 
150 Pa). The red line at 2 km indicates the current size of the nucleus, which is what we see today. For an initial radius of 2.2-3.0km, slightly 
larger than today, the pressure inside the nucleus exceeds the compressive strength of the layers we see today. 
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